Migration of bosonic particles across a Mott insulator to superfluid phase interface 
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We consider a boundary between a Mott insulator and a superfluid region of a Bose-Hubbard 
model at unit filling. Initially both regions are decoupled and cooled to their respective ground 
states. We show that, after switching on a small tunneling rate between both regions, all parti- 
cles of the Mott region migrate to the superfluid area. This migration takes place whenever the 
difference between the chemical potentials of both regions is less than the maximal energy of any 
eigenmode of the superfluid. We verify our results numerically with DMRG simulations and explain 
them analytically with a master equation approximation, finding good agreement between both ap- 
proaches. Finally we carry out a feasibility study for the observation of the effect in coupled arrays 
of micro-cavities and optical lattices. 
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Introduction - Effective many-particle systems in ar- 
tificial structures that can be well controlled in experi- 
ments have become an important testbed for the investi- 
gation of quantum many-particle and condensed matter 
physics. Extensive work with arrays of Josephsoni unc- 
tions [l| and ultra cold atoms in optical lattices Q has 
lead to substantial progress and seminal experiments. 

Very recently it has been shown, that arrays of cou- 
pled micro-cavities can host effective Bose-Hubbard mod- 
els for one Q and two [1] polariton components, related 
polariton models [B| and effective spin Hamiltonians 
The phase diagrams of these models @, HJ and the pos- 
sibility of a glassy phase have been discussed 0- As a 
new feature, this approach offers the possibility to con- 
trol and address individual lattice sites. Besides being a 
prerequisite for quantum information applications, this 
possibility opens the door to the study of many-particle 
systems which are inhomogeneous or out of equilibrium. 
In this work we exploit these strengths of local address- 
ability and controllability to study novel physical effects. 

We consider a one-dimensional Bose-Hubbard model 
of N sites, where sites 1 to Ni are in a Mott insulator 
and the rest in a superfluid regime. Initially there is 
no particle hopping between the two areas, which are 
both prepared in their respective ground states with on 
average one particle per site. At t = 0, we then switch 
on a small hopping rate between sites Ni and Nj + 1. 
As our results show, this causes all particles of the Mott 
region to migrate to the superfluid part, leaving the Mott 
part almost completely empty. 

We present time dependent DMRG simulations Q for 
finite systems with various parameters and an analyti- 
cal approximation using a master-equation for the case 
where the superfluid region is large and has no particle 
interactions. We find good agreement between both ap- 
proaches. These findings indicate that the scenario we 
consider could be used to experimentally study dissipa- 
tive quantum dynamics f]~oj ] where a part of the employed 



U » J 





part A 



part B 



t>0 



FIG. 1: Sketch of the considered scenario for a chain of 8 
sites with Ni = 4. Sites 1-4 form part A and are in a Mott 
insulator regime with U ^> J (symbolized by the narrow line 
that links the sites). Sites 5-8 form part B and are in a 
superfluid regime with U -C J (symbolized by the thick line 
that links the sites). At time t = 0, the hopping Ji between 
sites 4 and 5 is turned on. 



effective many body system acts as the heat bath, the 
properties of which can even be controlled and tested. 
Finally we discuss possibilities to observe the effect in 
both, coupled cavities 0] and optical lattices Q. 

Model and concept - We consider a Bose-Hubbard 
model, where the chemical potential, hopping rate and 
on-site interactions vary from site to site. The Hamilto- 
nian of our model with open boundary conditions reads, 



H 



N 



J j UX+i+h.c. (1) 



where ajj creates a particle in site j and rij = aijCLj- ftj, 
Jj and Uj are the chemical potential, hopping rate and 
on-site interaction at site j and Jn = 0. 

We consider a scenario (c.f. figure [1]) with Uj = U 
and fij = jjL for j = 1, . . . , Ni and denote this part of the 
chain part A. The remaining sites, which we will refer to 
as part B, have Uj = U and \ij = p, (j = Nj + 1, . . . , N). 
The hopping rates take values Jj = J for part A, i.e. for 
j = 1, . . . ,Ni — 1 and Jj = J for part B, i.e. for j = 
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FIG. 2: The particle density across the chain as function of 
time. Sites 1 to 20 start in a Mott insulator regime and the 
rest in a superfluid regime. U = 1.0, J = 0.1, /j, = 0, U = 0.2, 
(1 = 0, J = 1.0 and Jj = 0.1. 



TV/ + 1, . . . , TV — 1. while hopping between TV/ and TV/ + 1 
is initially zero, Jni = 0. Parts A and B are prepared in 
their respective ground states with filling factor 1, i.e. on 
average one particle per site in both parts. Hence, part 
A(B) initially contains TV/ (TV — TV/) particles. We assume 
that part A is operated in a Mott insulator regime, U ^ 
J, whereas part B is operated in a superfluid regime, 
U J. At t = 0, Jn, is then switched to a finite but 
small value Jat 7 = J/, where J/ <C U and J/ <C J. 

Numerics - We numerically calculated the initial 
state and the time evolution for chains of finite length us- 
ing the TEBD algorithm Q with matrices of dimension 



0.00 



- and a 4th order Trotter formula. 



20, timcstcps dt 
Truncation errors at each timcstcp were < 10~ 6 . 

As parts A and B are initially decoupled and both 
in their respective ground states, we calculate this initial 
state via an evolution in imaginary time, exp(-flH) with 
[3 — ► oo and H given by cq. ([T]) with Jjy 7 = 0. The imag- 
inary time evolution starts from a state with exactly one 
particle in each site, |1, 1, . . . , 1), and since it conserves 
the total number of particles in both parts independently, 
our initial state has unit filling with TV/ particles in part 
A and TV — TV/ particles in part B. 

We then simulate the time evolution, exp{—iHt), of 
this initial state, where H is given by eq. JT]) but now 
with Jjv, = J i- Figure [2] shows the evolution of the 
number densities in each site, (rij)(t) for a chain of length 
40. Sites 1-20 are in a Mott insulator regime with 
U = 1.0, J = 0.1, /i = and Jj = 0.1. For the remaining 
sites, U = 0.2, jl = and J = 1.0. The initial evolution 
is plotted separately with higher resolution in figure 
and shows that particles close to the boundary leave the 
Mott insulator region first. The total number of particles 
in part A, Na = Y^j=i( n j) 1S shown to decay to zero 



in figure 03. Initially TVa decays fast by one particle as 
the particle that is already next to the boundary leaves 
part A first. All further particles first need to travel to 
the boundary resulting in a slower decay of Na- Since 
part B has a finite size, one should expect recurrences in 
TV4. Whereas these appear on short time scales for small 
chains, the simulated time range in figure [5] is not long 
enough to see them. 

Analytical approach - Here we give an analyti- 
cal approximation for the case where part B has 
no interactions, U = 0, and is very large, N — 
Nj 1. To this end we derive a master equa- 
tion 11] for the reduced density matrix of part 
A, a = TrR(p), which reads, ^ = — i \H-i , a] — 



f* dsTr B {[Hi, [fl>(t - s), \GS)(GS\a}}}. Here, \GS) is 
the ground state of part B, Hi = — J)C^2r 1 ( a ] a J+ 1 + 
h.c), Hi = -Ji(a} Ni a Nl+ i + h.c.) and Hi(t) = 
with 



e iH t Hie -tH t with Hq = r //,(//, - 1, 

equation is valid up to second order in J and J/. We 
thus focus on the regime where U,J^> J,Ji- 

The Hamiltonian of part B is diagonalizcd via the 

/ 2 " y^jV-jVj • ( TTlj \ , 

Y N-Nt + 1 2^1=1 Sm \N-N I + l) ° l 

i( a ] a j+i + h- c = Y,iLi N ' Uib\bi 
with uji = — 2 J cos I Ar ^ , |. Then we find 



transformation a,j - 
so that — JY^,f=N 



N-N I + l / 

(GS\a Nl+1 (tW Nl+1 \GS) = (GS\a Nl+1 a^ Ni+1 (t)\GS) = 
(Jt) -1 J\ (2 J£), where Ji is a Bessel function of the first 



kind and we neglected terms of order 



These 



(N-Ni)* 

correlations decay sufficiently fast (cx t~ 3 / 2 ) and the 
time integral in the master equation can (after a trans- 
formation of the integration variable) be extended to the 
range (—00, 0) and the equation becomes an ordinary 



differential equation ^er 



-i{n Nl 



(nT n + mT m — in<d n + im<d r 
(r„+i + T m+ i - i6„+i + i6, 
jnjvj) is the state of site 



[Hi, a] 



\m Nl ) 



+ V»i + lvm+l 

1+1) f?i+l,m+l, an d 

TV/ with n particles and 
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a n . m — (nNr WImNi}- The energy shifts O, 
given by 6„ = ~{Jf/J)Xn for |x„| < 1 

e„ = 



are 
and 



— (J|/ J)Xtiv1 — Xri otherwise, whereas the 
decay rates T n are given by T n = Re f ( jf/ J) *J 1 — Xn J 

with x« = (U/J)(n -l)-(fj,- p)/(2J). There is thus a 
particle flow from part A to part B whenever \\ n \ < 1. 
One can identify two scenarios in which this flow is 
blocked: If the chemical potential is larger in part B, 
p > fx, all r„ are zero for fi — fl < —2 J. If on the other 
hand, the chemical potential is larger in part A, p < n, 
Ti is zero if fi — p > 2J, whereas other T n with n such 
that fi — p — 2 J > 2U(n — 1) remain nonzero. 

Hence, if part A is deep in the Mott insulator regime 
with U J and states |njyr) with njvj > 1 have very 
small occupation, the decay channels T n with n > 1 do 
not contribute and the particle flow vanishes whenever 
|/i — p\ > 2 J, that is if the difference of the chemical po- 
tentials in parts A and B is larger than the maximal en- 
ergy of any eigenmode of part B, i.e. — p\ > max/(|wj|) 
Nonetheless, for \i — p > 21 and moderate ratios U/J, 
there is a slow flow of particles from part A to B, which 
becomes more and more suppressed with increasing U /J. 
For U 3> J we thus approximate 0„ and T n by the values 
for n = 1 and obtain, 



da 

7/F 



-i Hi+ Qa^ Ni a Nl ,a 

+ T (2a Nl aa^ Ni - a^a^cr - era aw, ^ 



(2) 



We numerically tested the accuracy of eq. ((2]) for chains 
of N — 20 and N = 50 sites, where sites 1, ... ,4 form 
part A and the remaining sites part B. Figure 2] shows the 
results for various parameters. In the upper row, we plot 
the total number of particles in part A, Na = 
as given by the numerics whereas the lower row shows 
differences between Na as given by the numerics and Na 
as given by the master equation, [JV A ] numcriCB - [A^] master . 
We find good agreement between both approaches. 

The applicability of the master equation approach, 
which we have confirmed with our numerics, shows that 
the superfluid region behaves like bath and the Mott in- 
sulator like a quantum system under dissipation [l5| . Ex- 
perimental implementations of these scenarios would thus 
also allow to investigate processes as described by e.g. 
the spin boson model [l(| and its relation to the Kondo 
problem 17|. Moreover the properties of the bath could 
even be controlled and tested in these experiments. 

Experimental tests - As possible experiments we an- 
alyze here two realizations, polaritonic systems in cou- 
pled cavities and cold atoms in optical lattices. With 
the possibility to manipulate and measure individual lat- 
tice sites, coupled cavities suggest themselves as an ideal 
implementation of effective many-body systems for ob- 
serving the effect under consideration. One way to ini- 
tially separate part A and B is to generate a large, neg- 
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FIG. 4: Top row: The number of particles in part A, Na, as a 
function of time as given by the TEBD numerics. Bottom row: 
Corresponding differences between Na as given by the TEBD 
numerics and Na as given by equation ([2]), [ATA] numerics — 
[^master- U = 2, J = 1, U = 0.1 and N = 4 for all plots. 
The remaining parameters are J — 0.1, \i — jl = and TV = 20 
for plot A, J = 0.15, fj,-fl = and N = 20 for plot B, J = 0.1, 
jj, — p, = 1 and N = 50 for plot C and J = 0.1, fl — pi = 3 
and TV = 20 for plot D. For the parameters of plot C, good 
agreement between numerics and master equation is obtained 
for TV ~ 50 only because the density of states in part B is 
lower in the energy range of relevance here. 




FIG. 5: A Mott insulating and a superfluid region that are 
initially separated by a strongly negative chemical potential 
in site 4. U = 1, J = 0.1, U — 0.1, J = 1 and \ij = except 
for j — 10. Number densities {nj)(t) look similar as in figure 
H N A (t) (blue) and /no(i) (green). 



ative chemical potential in one cavity. This cavity thus 
stays empty and suppresses any particle tunneling or cou- 
pling between both parts. This chemical potential is then 
switched off and the particles start to migrate from the 
Mott insulator to the superfluid part. Figure [5] shows 
numerical results for such a scenario. Here sites 1 to 10 
are in a Mott and the rest in a superfluid regime. Site 10 
has initially a chemical potential /iio = —4 and thus no 
particle in it. As this chemical potential /xio is switched 
off during the evolution, particles start to flow into the 
superfluid part. The left plot shows number densities 
(rij)(t) whereas the right plot shows A^(i) (blue) and 
/iio(t) (green). In cavities the local number statistics can 
be measured using resonance fluorescence Q. 

A related experiment with cold atoms in optical lat- 
tices could be done by applying a magnetic field across 
the lattice. At first this field is constant and the system 
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FIG. 6: Particle migration resulting from a magnetic field 
gradient applied to an optical lattice with cold atoms, a: 
Uj{t), b: Jj(t), c: (n.j)(t) and d: S{k) at t = (blue), S{k) 



at t = 10 s (green) and S(k)\t=o - S(k)\t= 



(red). 



is prepared in the ground state of the corresponding ho- 
mogeneous model with Uj = Uq, Jj = Jq and fij = /io 
for all j. Then a magnetic field gradient is ramped up 
such that there is a spatially varying field across the lat- 
tice Since Uj/ 'Jj oc cxp(2y/Vj / E r ) (Vj is the lattice 
potential at site j and E r the recoil energy of the atoms) 
and furthermore Vj oc BJ 1 , where Bj is the magnetic 
field at site j, the ratio Uj/ Jj at each site can be tuned 
by the magnetic field. Hence a magnetic field gradient 
across the lattice can result in Uj/ Jj 3> 1 at one end of 
the lattice and Uj/Jj <C 1 at the opposite end. 

Figure [5] shows a numerical simulation for such an 



0, we set Uj 



Jj 



0.5 and 



experiment. At t 
fij = for all j (N = 20) and prepare the system in the 
corresponding unit filling ground state. Then we tune 
the interactions and hopping rates according to Uj(t) = 
0.5(l + 10- 4 (j-f )t) and Jj(t) = 0.5 (1- 10" 4 (j- f )t). 
The particles move towards the site with minimal Uj / Jj . 

Several quantities are accessible to measurements for 
cold atoms in optical lattices. Time-of-flight measure- 
ments Q reveal the number distribution in momentum 

space S(k) = iJ2f.i=i e27rkU ~ l)/N ( a ] a j)- Figured! 
shows the initial S(k) at t = (blue line), the final S(k) 
at t = 1000 (green line) and the difference between both, 
S(k)\ t= o — S(k)\ t= iooo (red line). Although the average 
ratio Uj/Jj remains constant, the peak broadens. In an 
inhomogeneous magnetic field, microwave pulses on field 
selective transitions can be used to selectively address 
atoms that experience a certain magnetic field and hence 
the density profile of the atoms can be measured 14 1 . 

Summary and outlook - We have considered an in- 
terface between a supcrfluid and a Mott insulator. Both 
parts were initially decoupled and cooled to their unit 



filling ground states. As a tunneling rate between both 
parts is switched on, all particles of the Mott region mi- 
grate into the supcrfluid part. This effect which we con- 
firm numerically by DMRG simulations and analytically 
with a master equation is caused by the high density of 
states of the superfluid region, whereas in Mott insulating 
area there is only one state in the accessible energy range. 
Our results also show how effective many body systems 
could be used to experimentally study dissipative quan- 
tum dynamics where the properties of the bath may be 
tested and cotrolled. Experimental observations appear 
feasible in arrays of micro-cavities and optical lattices. 
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